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Abstract 

The finite-size critical properties of the 0{n) vector ip"^ model, with long- 
range interaction decaying algebraically with the interparticle distance r like 
^-d-a ^ are investigated. The system is confined to a finite geometry subject 
to periodic boundary condition. Special attention is paid to the finite-size 
correction to the bulk susceptibility above the critical temperature T^.. We 
show that this correction has a power-law nature in the case of pure long-range 
interaction i.e. < cj < 2 and it turns out to be exponential in case of short- 
range interaction i.e. a = 2. The results are valid for arbitrary dimension d, 
between the lower ((i< = a) and the upper ((i> = 2a) critical dimensions. 

PACS numbers: 05.70.Jk, 64.60.Ak, 64.60.Fr 
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I. INTRODUCTION 



The critical behaviour at a second order phase transition depends upon the dimension- 
ahty (i, degrees of freedom ra, symmetry of the Hamiltonian (either in spin-space or in 
coordinate space) and interaction potentials. Generally speaking, the nature of the poten- 
tial of the model under consideration can describe different physical situations (for a review 
see reference The simplest interaction potential, which has attracted the attention of 
investigators, is the one corresponding to long-range ferromagnetic interaction decaying al- 
gebraically with the spin interdistance r as r~'^~'^, where d is the dimension of the system 
and a is the parameter controlling the range of the interaction. The interest in such type 
of interaction is tightly related to the exploration of the critical behaviour of systems with 
restricted dimensionality, in which no phase transitions occur otherwise. 

The investigation of systems with long-range interaction was initiated by Joyce in his 
paper on the ferromagnetic spherical model [Q. The results of Joyce were generalized to 
the 0{n) vector yj^ model by means of perturbation theory in combination with the renor- 
malization group techniques [^0] and the 1 / n-expansion |§] . These investigations were also 
extended to dynamical critical phenomena (see Ref. and references therein). Computer 
simulations also contributed in the exploration of the critical properties of such systems 
IPm^T^. The results of these simulations, obtained by means of the Monte Carlo method, 
concerned mainly systems with classical critical behaviour, in the sense that the critical 
exponents are those belonging to the Landau theory. Rigorous results were obtained for 
low dimensional systems with long-range interaction (see reference ITSlJT^ and references 
therein). 

The analytic exploration of the scaling properties of confined systems with long-range in- 
teraction took its starting point on the spherical model. The reason for choosing that model 
is the relatively simple nature of the mathematical expression entering the equations char- 
acterizing its thermodynamics (for a complete set of references on the subject see reference 
||T5|). Very recently these investigations were extended to 0{n) vector models [0,^, using 



the renormalization group approach and the e expansion to the one-loop order. In reference 
the Binder cumulant has been evaluated at the vicinity of the critical temperature. It 



has been found that the expression for that quantity can be deduced just by choosing an 
appropriate rescaling of the parameters in that evaluated for short-range interaction case. 
The authors of reference jl^ evaluated the susceptibility at the critical temperature Tc, as 
well as in the region, determined by the condition L/C, ^ 1, where L is the linear size of 
the system and ^ - the bulk correlation length. In this region the bulk critical behaviour 
dominates the finite-size critical one. It has been shown that the finite-size correction to the 
bulk critical properties of the system has a power-law nature. This result is distinct from 
that obtained for the case of short-range interaction, where the finite-size correction falls-off 
exponentially. 

In this paper we will investigate the finite-size scaling in a system with ferromagnetic 
long-range interaction at a fixed dimension d with d^ < d < d^. The parameters d^ = a and 
dy = 2cr are, respectively, the lower and upper critical dimensions of the model, with cr < 2. 
To this end we will use the approach developed in references [|T^-|20||. This method has been 
used successfully for the evaluation of critical exponents, as well as, critical amplitudes for 
various thermodynamic functions in the (p^ model with short-range potential. The most 
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important property of this approach is that the quadratic (temperature-dependent) term 
does not enter exphcitly the expansions, which can be used for both sides of the critical 
temperature T^. This method has been also used for the investigation of the theory of finite 



size-scaling ||21|-p3|. While a perfect agreement between the analytical results of reference 



2T| and those obtained by Monte Carlo simulations has been reported in reference ||24|| , 



references [0,^ showed disagreements with some of the known results. In particular we 
would like to emphasize their finding of the non-exponential decay of finite-size corrections 
to the bulk critical behaviour. 

We will consider spin system consisting of ra-component unit vectors associated with d- 
dimensional lattice and interacting via a pair potential of the general translationally invariant 
form. The Hamiltonian of this model is given by 



where is a short hand notation for the space dependent n-component field ^{x), tq = 
T^oc + ^0 (^0 T — Tc), a, b and uq are model constants. V = L'^ is the volume of the system. 
In equation (|1 . 1| ) , we assumed h = ks = 1 and the size scale is measured in units in which 
the velocity of excitations c = 1. We note that the operator is defined by its Fourier 
transform 



k 



The parameter (3 is set for the inverse temperature. Here we will consider periodic boundary 
conditions. This means 

V?(x) = L-"'^^(fc)exp(ifc ■ a;) , (1.2) 

k 

where fc is a discrete vector with components ki = 2Tini/L (nj = 0, ±1, ±2, ■ ■ ■ , i = 1, ■ ■ ■ ,d) 
and a cutoff A ~ {a is the lattice spacing). In this paper, we are interested in the 
continuum limit i.e. a — > 0. As long as the system is finite we have to take into account the 
following assumptions L/a— >oo,^— s>oo while C,/L is finite. 

The critical behaviour of the model Hamiltonian ( p..l|) , has been investigated in details 
in the early 70 's. The main focus of interest has been turned to the evaluation of the critical 
exponents. Using renormalization group techniques, it has been shown that the critical 
behaviour of this model is dominated by the long-range interaction for < a < 2 l^,^. In 
this case the critical exponents are a dependent, in particular the Fisher exponent r^^ = 2 — a. 
As long as a becomes of the same order as 2 — 772, where t]2 is the Fisher exponent of the 
short-range model, a crossover from the long-range critical behaviour to the short-range 
critical one takes place @,PJ3l- ■^o'^ > 2 — r]2 the critical behaviour is dominated by the 
short-range interaction and the critical exponents are those of the pure short-range model. 
Another way to establish the relevance of the long range term for a < 2 — 772 is presented in 
reference 



The plan of this paper is as follows. In Section || we describe the renormalization scheme 
for the bulk theory with long-range interaction. In Section |ITT| we discuss the effects of 
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confined geometries on the bulk critical behaviour. We investigate the influence of the long- 



range interaction on the flnite-size correction to the bulk critical behaviour. In Section ^ 
we discuss our results briefly. An appendix is added in order to complement the results of 



Section III 



II. THE BULK SYSTEM 
A. Bare theory 

For simplicity here we will consider the model with pure long-range ferromagnetic in- 
teraction i.e. the parameter a controlling the range of the interaction is smaller than 2, 
value characterizing the short range interaction. In other words, we will consider the model 
( |1.1| ) with the parameters a = and 6=1. We believe that one has flrst to understand 
this model before starting to explore the model with both long and short-range interactions 
present. The investigation of the critical properties of the complex structure of this model 
will presented elsewhere. 

It has been shown that the Hamiltonian with pure long-range interaction can be treated 
with fleld theoretical renormalization group techniques |^. The renormalization constants 
as well as the fleld theoretic functions were calculated. A superflcial discontinuity of the 
anomalous dimension of fleld theories occurs as soon as o" = 2. This however is true only at 
this particular point, and as long as we are far from that point, we can use the model with 
the spectrum ro + In this case the renormalization constants, to one loop order, are 
given by: 

Z^ = l + 0{u^) (2.1a) 
n + 2 

Zr = l + M + 0{u^) (2.1b) 
= l + ^^M + 0(u2). (2.1c) 

Here, as usual, is the scaling fleld amplitude, Z^ the coupling constant renormalization, 
and Zr the renormalization of the y?^ insertions in the critical theory. The parameter e = 
2a — d denotes the deviation from the upper critical dimension. 

The application of the method proposed by the authors of references ||TH|-^ to systems 
with pure long-range interaction can be established easily following the way this has been 
done in combination with the e expansion ||^. In particular we turn our attention to the 
inverse bare susceptibility, related to the two-point vertex function, at flnite external wave- 
number k 

-^ = r(2)(fc,ro,no,A,rf). (2.2) 

This function is, of course, well deflned at T ^ for dimensions d > a. The parameter 
To is, as usual, a linear function of the reduced temperature t = {T — Tc)/Tc. Its critical 
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value is determined by the condition x(0) ^ = O5 from which we obtain a natural implicit 
definition for the function 

roc = roc{uo,A,d). (2.3) 

Following reference [0, we will consider the difference tq — tqc, which is a function of the 

correlation length (i.e. vq — tqc = h{^, uq, d)), instead of tq itself into the expressions for the 



different vertex functions. So, we will consider the vertex functions F*^^) as depending on 

the parameters tq — tqc, Uq, A and d. We will denote the dimensionally regularized vertex 


function by F''^^ = F^^-* (ro — roc, mq, c^)- Correspondingly the dimensionally regularized 
critical parameter tqc will be denoted by rQc{uo,d). From simple dimensional arguments, 
one obtains the relation 

roc oc Mq^^ (2.4) 

The method used here differs from the ones when the e expansion come into play, by the fact 
that the critical parameter roc is non-vanishing within the dimensional regularized theory. 
It is clear that the use of the e expansion implies an apparent vanishing of tqc, as seen 
formally from the relation (|2.4|) between tqc and uq, which for infinitesimal e does not yield 
a contribution at finite order to the perturbation theory. 

In order to make clear the definitions introduced along this section, we will present here 
the one-loop results. For the two point vertex function we get 

-s/a- 

F(2) =ro + k''-{n + 2)uoAd/-^ + 0{ul), (2.5) 
where the geometrical factor Ad^a is defined by 

The four point vertex function is 

~e/a 

g y 
a J e 



?W =u,-{n + %)ulAd,, (1 - i) !l + o{ul). (2.7) 



In the next subsection we will discuss the renormalization of the bare theory for fixed 
dimension d confined between the lower and the upper critical dimensions given by a and 
2(7, respectively. 



B. Renormalized theory 

It is well known, that the perturbative results of the bare theory do not provide a correct 
description in the critical region ^ — > cx), for dimensions below the upper critical dimension. 
This can be outwitted by taking advantage of the ideas of the renormalized theory, which 
furnishes a mapping from the critical region to non-critical one. Our starting point is the 



5 





expression for the A^-point vertex functions Tq {^,Uo,d) obtained from T^'^^ (ro — roc,Uo,d) 
by switching from the variable tq — tqc to ^. This is possible because of the fact that 
the reduced temperature is tightly related to the correlation length. The deviation of the 
parameter tq from its critical value will be a fixed quantity and hence the correlation length 
will have the same property. 

We treat the theory by using the minimal subtraction scheme at fixed dimension d. To 
this end, we introduce the renormalized quantities 

ifix) = Zy\n{x) (2.8a) 
'^o = '"oc + ZrT, ro - roc > (2.8b) 
^d,aUQ = n'ZuZ~'^u. (2.8c) 

Then the renormalized vertex functions takes the form 

fW(e,«,/i,rf) = Z^/'T^i^,i^^Z.^Z^\A,^^,d). (2.9) 

In addition, we require that the renormalization constants and Z^ absorb just the poles 
of at the upper critical dimension, which turns out to be 2a for the model under consid- 
eration. In equations ( |2.8| ) and ( p.9|) the parameter fi is an inverse reference length, which 
will be chosen as the amplitude of the asymptotic bulk correlation length. 

In the following we proceed in a standard way, by deriving a differential renormalization- 
group equation for the vertex function r*^^). This is achieved by taking the derivative of 
equation ( |2.9| ) with respect to fi at fixed constant Uq and tq — roc- This leads to 

+ Pu{u, e)du + ^iVC^(M)]f(^)(e, u, fx, d) = 0, (2.10) 

with 

Pu{u,e) = {fiduu)o, (^{u) = {fidf,\iaZ;^^)Q. (2.11) 

Here the subscript indicates that the differentiation is performed at fixed parameters of the 
bare theory. Using the method of characteristics, a formal solution of the renormalization- 
group differential equation (|2.10|) is given by 



f(^)(e,«,/x,d) = f(^)(e,«(£),£/i,ci)exp(^-y^ UnY) ■ (2.12) 

Here ({i){u{£),d) and u{i) is the solution of flow equation 

i^ = PuHi),d\ (2.13) 
with the initial condition u{i) = u. The most convenient choice for the flow parameter I is 

^^i = C'■ (2.14) 
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The renormalized vertex functions can be written as 



fW(e,«,/i,^i) = r'+^'^-'^^^/^^H/ie,«,c^), (2.15) 

where the amphtude functions f^^^ are dimensionless. The renormahzabihty of the (p^ model 
with long-range interaction for dimensions d less than the upper critical dimension 2cr is a 
warranty for the finiteness of the F^^-* at fixed u, jj, and ^. From ( p.l2[ ), ( |2.14| ) and ( p.l5D , 
we obtain the expression 

u, d) = d) exp [ - / Un^) . (2.16) 



As usual, equation (|2.16| ) lies in the basis of the mapping of the amplitude function 
f^'^\n^,u,d) from the critical region, where the perturbation theory breaks down to the 
noncritical region, where the perturbation theory is applicable. 

C. Asymptotic regime 

In the asymptotic limit, determined by (£ — 0, ^ — oo), the coupling u{i) approaches 
the fixed point u* = u{0), which is the zero of the f3 function of the theory i.e. 

Puiu*,0,d)=0. (2.17) 

For ^ — > oo, equation ( p.l6|) , takes the asymptotic form 



/(^)(/ie,«,rf) ~ AW/(^)(l,«^o,rf)(/iO^''/^ (2.18) 

with the critical exponent t] = —Cip{u*, 0, d) and the nonuniversal amplitude 

^^^^ = exp I'lUn - UO)]y) ■ (2.19) 

The results obtained here can be applied very easily to the particular case of the bare 
susceptibility corresponding to the particular value N = 2. This is given by 

X = Z^iu, d)e[f^'\l, uii), d)]-' exp (^1' Un^) ■ (2.20) 
The critical behaviour above Tc reads 

X = Ap~\ (2.21) 

where 

= eoZ^iu,d)[A^'^f^'\l,u*,0,d)]-\ (2.22) 
Here we have used fx = and the asymptotic form ^ = " with 7 = i^(2 — r]). 
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III. THE FINITE SYSTEM 



A. The effective Hamiltonian 



For the evaluation of the effective Hamiltonian of the finite system, we need to calculate 
the free energy per unit volume. This is defined by 



f = -L~'^lnZ; Z= Vipexp{-pn). 



Here we are also interested in the expression of the susceptibility, which is defined by 

1 



Z 



(3.1) 



(3.2) 



Following references p6|,p7|, we split the field = $ + S into a mode independent part 



(3.3) 



which is equivalent to the magnetization and a part depending upon the non zero modes 



E = L-'^^V(fc)exp(ifc-a;). 



(3.4) 



Consequently (|TTT|) is decomposed as 
with the zero-mode Hamiltonian 



Whence the partition function takes the form 

poo 

Z= / d^expl-iHo + m^m 

J — oo 

where 



r($^ 



In / VEexpl-Hi] 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



contains the contribution from higher order. 

Now instead of the decomposition of reference we will use the modified perturbation 
theory proposed in reference |^T|]. There, an appropriate decomposition of the 0{n) vector 

model with short-range interaction has been presented. This way has been proven to give 
a good quantitative results above, as well as below the critical temperature Tc. Applying 
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that method to our model we find, that the higher-mode dependent part Hi can be spht 
into 

^1 = 11 d'4iro + {n + 2)u,Ml)T? + [V'^'^T.f] 

+ ^ j d''x[?,u,{^^ - Ml)T? + 2mo$S3 + l^oS^], (3.9) 

where we have introduced the magnetization 

M2 = (02). (3.10) 



Further a diagrammatic expansion of r($^) can be represented by two point vertices propor- 
tional to Mo($^ — Mq) in addition to the three and four point vertices ~ no$ and ~ uq. In 



this way, the finite-size perturbation theory is obtained as an expansion of r($^) in powers 



of Mo($^ — Mq), Mo and Mq^^. In the following, instead of r($^) we consider L~'^r($^) which 
remains finite in the limit L — oo. To the leading order the finite-size correction reads 



k 

-^^uli^' + 2Mie)S,{T,) + OK, ulct?, ul{^' - Mlf), (3.11) 



where 

Sm{r) = L-''Y^{r + {n + 2)u^MI + fc")"'". (3.12) 

k 

Here we will investigate only the high temperature regime i.e. T ^ T^, so we will set 
Mq = 0. a rearrangement of equations (|3.6|) and ( p.ll|) leads to the effective Hamiltonian 

KfT(ro, Mo, L, $2) = + + 0($6)]^ (3.13) 

where 

i? = ro + (n + 2)uQL-'^Y^{rQ + fc")-\ (3.14a) 

k 
k 

Finally, the free energy ( p.l|) reads 

/ = L-4(0) - L-*^ In Zeff, (3.15a) 

poo 

Zefr= / rf<l>exp[Kff(ro,uo,L,$2)] (3.15b) 
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and averages such as are calculated from ( p.2|) as 

1 r°° 

<$2p^= / rf$$'^exp(-Kff(ro,Mo,L,$2)). (3.16) 

Let us note that for the purpose of comparison with numerical results it is convenient to 
keep the exponential structure of the integrand in the various thermodynamic functions. 
Notice that an appropriate rescaling of the field $ in equation ( p.l6| ), permits to write the 
momenta ClS db function the 'scaling variable' 

z = RL'''^U-^'\ (3.17) 

Whence 

/$2P\ = (m^Z'^)~p/2 Jo "-^-^ exp ^zx ) 

dxx'^~^ exp (^~zx'^ — \x'^) 

These functions can be expressed in terms of the confluent hypergeometric function (see 
Appendix 0). 



B. Renormalization 



As it has been mentioned in the previous section, in this paper, we will use the approach 
of references [|18|- pO|] . Here we will consider the limit of an inflnite cutoff i.e. A — oo at fixed 



~ "^Qo where roc is the bulk critical value of tq. Let us recall that the advantage of using 
this method is its direct application to fixed space dimensionality without resorting to 
the e expansion. For systems with short-range interaction this method has been tested and 



accurate results, for various thermodynamic functions has been, obtained |^T|. Since roc is of 
order 0{u^^'^) in dimensionally regularized form, we find it safe to replace in our expressions 
tq by the difference tq—tqc- The effective Hamiltonian will be function of tq—tqc but not of ro 
itself. The effective constants R and U defined in ( p.l4| ) are also functions of this deviation, 
which itself is a function of the correlation length. The perturbation approach consists of an 
expansion of with respect to the renormalized counterpart u of the coupling constant Uq 
at fixed dimension. Since the finite nature of the geometry of the system does not alter the 
ultraviolet divergences, the usual bulk renormalizations are enough to describe the scaling 
properties the finite system. 

Here we are interested in the evaluation of the explicit expressions of various thermody- 
namic quantities for the finite system. Different crossover regions, in terms of renormalized 
parameters, are under consideration. The corresponding renormalized effective coupling 
constants to equations ( |3.14| ) take the form 
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rji) 



-e/cr 



A. 



(3.19a) 



-e/cr 



^ 1 ^ r(£) 



-e/cr 



where we have used the definitions [|T7| 



Siir) = -Ad,.e-\'-'/^ + L'^-'^FdArLn. 

S2{r) = Ad,^e-' (l - ^) r-^/'^ - L'^-'^^F^ArLn, 



(3.19b) 

(3.20a) 
(3.20b) 



and 



Fd,a {y) 



1 



{2'kY 
Here the function 



dxx 2 ^ii^. 



(7 (7 

2 ' 2 



ct/2 



{2'nY 



(3.21) 



k=0 



r {ak + 13) 



(3.22) 



is the so called Mittag-Leffler type functions. For a more recent review on these functions 
and others related to them, and their application in statistical and continuum mechanics 
see reference . A brief summary of some of the properties of the Mittag-Leffler functions 



are presented in reference [17 



By making the choice such that the flow parameter £ = £(t, L) satisfies the following 
relation 



in the critical region, equations ( 3.19|) imply the scaling forms 



(3.23) 



(3.24) 



Using the asymptotic expressions for R{i) and U{£), we obtain the asymptotic form of the 
variable z, defined in ( p.l7| ). It is given by 



-1/2 



(3.25) 



From this expression, one can convince himself, easily, that the function z{tL^^'^) has an 
expansion in terms of \/u*. 
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C. The Susceptibihty 



In this section, we present our resuh for the finite-size scahng function of the magnetic 
susceptibihty. Here we give the general expression, obtained by expanding with respect to 
the coupling constant u about its fixed point u*. The dimension is kept fixed i.e. we are not 
going to expand in the vicinity of the upper critical dimension as it has been done earlier 
||17|| . Furthermore, we won't expand the exponential weight of the integrand in equation 
( p.l6| ). In this way consistency with the correct one- loop bulk expressions is ensured in the 
bulk limit, reached by sending the size L of the system to infinity. As it has been mentioned 
in the previous sections the arbitrary reference length of the renormalized theory will 
be chosen as the amplitude of the correlation length i.e. fi^^ = ^q. 

In reference it has been shown that, for periodic boundary conditions, the ultraviolet 
divergences of the susceptibility x of finite systems in the continuum limit i.e. a — >■ are 
identical to those of the bulk susceptibility Xb- In reference [O], the application of this 



method has been extended to systems with long-range interaction and some results for the 
susceptibility in the asymptotic regimes L/C, S> 1, as well as L/^ ^ 1 has been obtained. 
The results obtained there were limited to dimensions close to the upper critical one. Here 
we will extend these results to the whole interval of dimensions between the lower critical 
dimension d^= a and the upper critical one dy = 2a. 

The average (0^), entering the definition of the susceptibility ( p.2|) , is defined with the 



statistical weight exp[— (7-^o + T)] given in equation ( p.7|) . Obviously, in the bulk limit we 
recover the bulk susceptibility i.e. limi_,oo X = Xb- 

For systems confined to a finite geometry, the renormalized susceptibility xr^ ^is a func- 
tion of ro — roc , uq, and L, can be introduced as 

L, /i, d) = Z-'^xi^, ifZuZpAd^^u, L, d), (3.26) 

where Z^, and Zu are the bulk Z amplitude factors defined in Section |T|. 

A renormalization-group equation for xr is obtained by deriving expression ( p.26|) for the 



susceptibility with respect to the parameter /i at fixed tq — tqc (function of the correlation 
length), uq and d. Since the linear size L of the system does not renormalize, we get 

[fidf, + (3udu - Q Xr{^, u, L, /i, d) = 0, (3.27) 

where the renormalization group functions l3{u) and ({u) are defined in equation ( p.ll| ) for 
the bulk theory. A formal solution for this equation is given by 

Xr{C, u, L, /i, d) = xr{^, u{i), L, ifi, d) exp (^^ U^')'^^ , (3-28) 

where the parameter i can be chosen arbitrarily. The most convenient choice is that, for 
which i satisfies the relation ( p.23|) . Remark that the renormalization group equations (p.lO|) 



for the bulk system and (|3.27|) for the system confined to the finite geometry are similar. 
However because of the finiteness of the size L of the system a careful consideration of equa- 
tions (|3.27| ) and its solution (|3.28|) is in order. In this case, we introduce the dimensionless 



amplitude function f^iz) according to 

XRiL, ^, u, /i, d) = L'^f^ifiL, ij,^, u, d). (3.29) 
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In the asymptotic regime, given by ^ <^ 1, we obtain from the formal solution (|3.28| ) of the 
renormalization group equation 



Xr{L, e, u, /X, d) ~ LH^ u\ d). (3.30) 

Using the fact that and are functions only of the ratio L/^, which follows from 
equation ( ^.23] ), we can write the bare susceptibility x = Z^XR: as well as, the renormalized 
susceptibility xr iii the following scaling form 

Xfi(L,e,u,/x,(i)~L2-v' 



1 + (lyP r (I) , (3.31) 



where Y{z) is a scaling function of its argument. 

Equation ( p.31j ) is the final expression for the susceptibility of the finite system. This is 
the complete expression in whole — plane. In the following we will consider different 
regimes with respect to the ratio C,/L. For this purpose we investigate the 0{n) vector (p^ 
model. Here we turn a special attention to the limiting case ^/L ^ 1. 

In the remainder of this section we will investigate the behaviour of the susceptibility 
using an ordinary perturbation theory. Accordingly, the standard one-loop expression for 
the inverse susceptibility above the critical point for a system confined to a finite geometry 
reads 

= ro + (n + 2)uoL-'' J^i^o + fc")"' + 0{ul). (3.32) 

k 

The susceptibility for the bulk system is obtained by sending the size of the system to infinity 
i.e. the sum in the last equation is replaced by integrals. It is 

Xb' = ro + in + 2)uo J d^k^r^ + k'^Y^ + 0{ul). (3.33) 

After some straightforward calculation one gets 

X(t, Mo, L, d)-' =t[l + {n + 2)uoL^''-\tL'')-^ (1 + tL'' F^^^itL''))] . (3.34) 

From this expression, we see that it is not allowed to set y = 0. In other words to take the 
limit t — > 0, while the size of the system is kept fixed. In the opposite limit L/^ ^ 1, the 
right hand side of equation ( |3.34| ) is well defined and the result gives the finite-size correction 
to the expression of the bulk susceptibility. However, the behaviour of the function F^^ain) 
is strongly dependent upon the value of the parameter a. Indeed, it depends upon the 
nature of the interparticle interaction in the system. Here we will discuss the influence of 
the interaction on the behaviour of the finite-size correction. 



The function Fii^„{y) has the following large y asymptotic behaviour |T7 
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for the case < cr < 2, and 



y 



(3.36) 



for the particular case a = 2. These results show that the last term in equation ( |A9| ) is just 
cancelled by the first term in equations (|3.35|) and ( |3.36|) . 

In the case of long-range interaction < a < 2, we obtain for the susceptibility, after 
renormalizing the theory, 



X = Xb 



1 - u*{n + 2)2'^7r-'^/2 



-d—a 



+ 0{u* 



(3.37) 



in agreement with the finite-size scaling hypothesis. Equation (p.37|) shows that the finite- 
size scaling behaviour of the system is dominated by the bulk critical behaviour, with small 
correction in powers of L. First the power law fall-off of the finite-size corrections to the bulk 
critical behaviour, due to long-range nature of the interaction, was found in the framework 
of the spherical model |]30| , |3T| , which is believed to belong to the same class of universality 
as the 0{n) vector model with n oo. 

It should be noted that the above result ( |3.37| ) cannot be continued smoothly to the case 
of short-range interaction cr = 2. In this particular caseFrf_2(l/) (see equation (p.36|) ) falls 
off exponentially fast and, correspondingly, the finite-size corrections to x ci-re exponentially 
small: 



X = Xb 



l-u*{n + 2)d{27itL 



2n(1 



^d)/4^-LVt 



(3.38) 



23[| . Notice that by expanding to the first order 



This result was obtained first in reference 

in e = 2cr — c? in equation ( ^.37] ) and ( pj.38D , we recover, the corresponding results (3.10) and 



(3.31) of reference 



Here an important remark is in order. In appendix we show that the results ob- 
tained for the susceptibility in this section are identical to those obtained, using the mode 
expansion. This is in disagreement with the conclusions of reference claiming that the 
mode expansion is inadequate for the description of the finite-size scaling in the region of the 
phase diagram determined by the condition L/C, ^ 1. Especially for the case of short-range 
interaction this result is in agreement with that obtained by means of Monte Carlo. 



IV. DISCUSSION 



In this paper, we have investigated the finite-size scaling properties in the 0{n)- 
symmetric ip'^ model with long-range interaction potential decaying algebraically with the 
interparticle distance. By means of the method, which consists of the use of the minimal 



subtraction scheme applied to a fixed space dimensionality developed in references [p!^-|20 



for the bulk system and in [^I|-|23[| for the one with finite linear size, we generalized the 
results of reference [l7| obtained for the susceptibility by means of the e = 2a — d expansion. 
The obtained results are in a full agreement with those obtained in reference |jl7 
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Here we restricted our calculation to the critical domain T ^ Tc and investigated the 
behaviour of the susceptibihty to the one-loop order in the coupling constant Uq of the bare 
theory. We have turned our attention to the special case, in the phase diagram, where the 
bulk critical behaviour is dominating the finite-size scaling one. We have found that the 
finite-size correction falls off algebraically ( p.37| ) in the case when a < 2 and exponentially 
( p.38| ) in the particular case cr = 2, characterizing the short-range interaction. These results 
are obtained by using two different methods: (i) The standard mode expansion of references 
26| , |27[| and (ii) The ordinary perturbation theory used in references p2| , p3| . We would like 
to mention here that our results does not agree with those of reference [E3| claiming that the 
method built on the mode expansion does not adequately describe the critical behaviour of 
the finite system. 

Note that by evaluating the Binder's cumulant ratio, we obtain a result, which have the 
same behaviour as a function of e as that of the result of references IITgIJTTII. It seems that 



to the first loop order, the method used here does not ameliorate the result obtained by 
means of the e-expansion. It is possible that higher loop order can improve the result in 
comparison to the Monte Carlo method of reference |]16| . 

In this paper, we concentrated our attention to a field theoretical model in the continuum 
(scaling) limit. Consequently the cutoff is send to infinity. By accounting the finite cutoff 
effects we expect that the finite-size scaling will be violated in a similar way to the case 
short-range interaction p2|. Note that the same effects can be obtained if one takes the 



model ( |1.1| ) with the parameter a controlling the range of the interaction to be larger than 
the value, a = 2, characterizing the short potential. This would verify whether the results 
obtained in references at the spherical limit remains true for finite n. 

Another possible extension of the results obtained here, in the static limit, is the appli- 
cation of the finite-size scaling theory to systems including critical dynamics. This can have 
direct implication to so systems exhibiting quantum critical behaviour. 
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APPENDIX A: FINITE-SIZE CORRECTIONS TO THE BULK SUSCEPTIBILITY 



In the appendix, we evaluate the finite-size correction to the bulk susceptibility in the 
region L/^ ^ 1 for arbitrary dimension of the system. Here we use the mode expansion and 



we will work in the one loop order as explained in Section |T|. We are interested in particular 
in dimensions d, such that a < d < 2a. In the leading order of the non-zero (fc ^ 0) modes 



the effective Hamiltonian ( 3.13 ) reads 



'Hcs{ro,uo,L, 



(Al) 



where 
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R = rQ-roc + {n + 2)uo 



k 



U = uo-{n + 8)ulL-^Y^{ro - r,, + k^^Y 



(A2a) 
(A2b) 



We have incorporated here the finite shift roc = —{n + 2)uq jj^ k '^ + 0{uq) of the parameter 
tq at one loop-order. Let us recall that in general one does a double expansion: one in 
the modes and the other in the coupling constant mq. In other words the fact that we are 
working in fixed space dimension does not mean that the parameter uq should be kept fixed 
or we are not allowed to expand with respect to the small parameter uq. This observation 
will be of a great importance in the following. 

In the one-loop order the difference tq — tqc is proportional to the reduced temperature 
t i.e. To — roc = at. In the following we will choose the coefficient a = 1. In terms of t, 
equations (|A2|) reads 



R 

U 



t- (n + 2)MoAi(t) 

,2 



where 



uq- {n + 



u, 



/ (t + k 

k 



cr\-2 



+ {n + S)ulA2{t), 



k 



kn 



(A3a) 
(A3b) 

(A4) 



Now we will evaluate the susceptibility for the system confined to a finite geometry. In 
the present approximation it is 



X 




with 



dxx"'~^^ exp 



1 



(ixx" ^ exp (- 



(A5) 



(A6) 



The integrals appearing in equation ([A6| ) can be expressed in terms of the confluent 
hypergeometric function t/(a, 6; z) according to 



dyy 



2 





1 












2' 





for 2; > 0. 



(AT) 



This function has a well known analytic properties (see reference ||33|| ). 

In the region tL^^" ^ 1, corresponding to 2; ^ 1, using the asymptotic form of the 
function Q{z) for large argument, which follows from that of the confluent Hypergeometric 
function, we obtain for the susceptibility 



X 



1 



1- (n + 2) 



U 



(A8) 



16 



Here we would like to mention that this result is not the final expression for the suscepti- 
bility of the finite system. Here it easy to make a mistake by taking, to the lowest mode 
approximation, the coupling constants R and U to be equal to the initial coupling constants 
^0 ~ foc and Uq, respectively. One has to keep in mind that apart from the lowest mode 
approximation, a loop expansion in mq is involved in the calculations. So, here one must 
take into account the full expressions of the constants R and U up to the lowest (one-loop) 
order in uq. Consequently by expanding in equation ( [X^ ) to order 0{uq) we arrive to the 
final result for the susceptibility 



x = t~ 



y 



(A9) 



Here we have used y = tL'^ and Ai(t) = —L'^"~'^Fd^fj{tL''). As one sees equations ( |3.34| ) and 
( |A9| ) are similar and all what it has been written after equation (|3.34]) remains true here. 
The conclusion we can draw from here is that both the ordinary perturbation theory and the 
method based on the mode expansion give only one and the same result for the susceptibility 
in the region, where the bulk properties of the system are dominating the critical behaviour 
of the system i.e. in the region L/^ ^ 1. 
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